Radiation-Induced Magnetoresistance Oscillations in a 2D Electron Gas 
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Recent measurements of a 2D electron gas subjected to microwave radiation reveal a magnetoresis- 
tance with an oscillatory dependence on the ratio of radiation frequency to cyclotron frequency. We 
perform a diagrammatic calculation and find radiation-induced resistivity oscillations with the cor- 
rect period and phase. Results are explained via a simple picture of current induced by photo-excited 
disorder-scattered electrons. The oscillations increase with radiation intensity, easily exceeding the 
dark resistivity and resulting in negative-resistivity minima. At high intensity, we identify additional 
features, likely due to multi-photon processes, which have yet to be observed experimentally. 

PACS numbers: 73.40.-c, 78.67.-n, 73.43.-f 
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The electrical transport properties of a 2D electron 
gas (2DEG) in a perpendicular magnetic field have been 
studied extensively over the past two decades in con- 
nection with the quantum Hall effects. However, re- 
cent experiments, in which such systems are subjected to 
electromagnetic radiation, reveal a surprising new phe- 
nomenon. The initial experiments of Zudov et al. 0, 
as well as more detailed subsequent studies conducted 
by Mani et al. 2} and Zudov et al. show that a 
peculiar oscillation of the longitudinal resistance is in- 
duced by the presence of microwave radiation in systems 
at high filling factor. Unlike the familiar Shubnikov-de 
Haas (SdH) oscillations which are controlled by the ratio 
of the chemical potential, fi, to the cyclotron frequency, 
ui c , these radiation- induced oscillations are controlled by 
the ratio of the radiation frequency, u>, to the cyclotron 
frequency. According to Ref. |2], the minimum resistance 
values are obtained near oj/uo c — integer + 1/4. These 
results are surprising since naively one would only ex- 
pect a peak at u> — u> c due to heating at the cyclotron 
resonance. Interest in this phenomenon was heightened 
when, using high mobility samples, both Mani et al. Q 
and Zudov et al. observed that as the radiation inten- 
sity is increased, the minimum resistance values approach 
zero and give rise to new zero-resistance states. In this 
Letter, we present calculations which reveal the origin 
of the resistance oscillations, explain their period and 
phase, hint at the manner in which zero-resistance states 
may develop, and predict higher-order effects, likely due 
to multi-photon processes, which set in at high intensity. 

It is instructive to begin by exploring a crude treat- 
ment through which a simple physical picture can be es- 
tablished. The electronic states of 2D electrons without 
disorder in a perpendicular magnetic field, B, are the 
Landau levels, e n = nuj c , where lo c = eB/m*c, the n 
are nonnegative integers, and m* is the effective mass. 
(Here, and always, we set h = 1 and drop the constant 
energy shift of uj c /2.) Each level has a degeneracy per 
unit area of l/2ir£ 2 where £ — y c/eB is the magnetic 
length. In the presence of disorder, the Landau levels 
are broadened. In correspondence with experiment, we 




FIG. 1: Simple picture of radiation-induced disorder-assisted 
current. Landau levels are tilted by the applied dc bias. Elec- 
trons absorb photons and are excited by energy ui. Photo- 
excited electrons are scattered by disorder and kicked to the 
right or to the left by a distance ±Ax. If the final density 
of states to the left exceeds that to the right, dc current is 
enhanced. If vice versa, dc current is diminished. Note that 
electrons initially near the center of a Landau level (where the 
initial density of states is greatest) will tend to flow uphill for 
uj/u! c ~ integer + 1/4. 



consider states occupied up to high Landau levels (the 
filling factor is about 50 for the experiments of Ref. 

EJ. 

We suppose that the normalized eigenstates in the pres- 
ence of disorder are ip a (r), with eigenvalues e Q . Now 
we turn on the microwave radiation. When an electron 
absorbs a photon, it is excited by an energy to. In the 
absence of disorder, conductivity is independent of the 
distribution of electrons over Landau levels (Kohn's the- 
orem '4]) so photo-excited electrons can make no addi- 
tional contribution to the dc current. In the presence of 
disorder, these electrons can be scattered by impurities 
as they are excited. We use a generalization of Fermi's 
golden rule to obtain a position-dependent rate. To sec- 
ond order in the radiation, the rate at which an clcc- 
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tron is scattered from initial state a' to final position r 
is w a i(r) = 2nJ2 a \^a(r)\ 2 5(e a ~ - w) | (a|e^x|a') | 2 , 
where (a|eSx|a') is the matrix element of the electric 
dipole operator of the radiation field. Initially, the 
eigenstates are occupied with Fermi distribution func- 
tion np(e a >) about chemical potential /i, and we insert 
the factor [n.i?(e Q /) — jif^c*)] so that final states must 
be unoccupied. Then the additional current density 
due to the radiation is AJ X (R) = — e J decPAr [np(e) — 
n F (e + u)] £ a , IVV {r')\ 2 S{e - e a >)w a > (r) Ax. Here R = 
i(r + r'), Ar = r — r' = (Ax, Ay). For the average over 
disorder (denoted by an overbar throughout this paper) , 
we write this expression as 

AJ X (R) = -2ire J ded 2 Ar [n F (e) -n F (e + u)] 

xN(r,e + u)N(r',e)M(r,r')Ax (1) 



where N(r, e) = J2 a IV"a(r)| 2 8(e — e a ) is the average lo- 
cal density of states, which in high Landau levels obeys 
N(e + nuo c ) = N(e). M (r, r') > is defined by extract- 
ing these density of states factors, and is assumed to be 
independent of the energies e, e + w, and to be a func- 
tion of |Ar| only. To measure the conductivity in linear 
response, a small dc electric field, Ed c , is applied in the 
x-direction. This tilts the energy spectrum as depicted in 
Fig. ^ and thus we assume that the only effect on Eq. (|TJ 
is to make N(r, e) = N(e — eEd c x). AJ X vanishes to zc- 
roth order in Ed c by symmetry, and for linear response, 
we expand to first order in Ed c , and divide by Ed c to 
find the longitudinal conductivity, a xx . We find that the 
radiation-induced change in the longitudinal conductiv- 
ity is proportional to an integral of the partial derivative 
(d(N(r, e+u)Af{r', e))/dAx) R . The density of states can 
be roughly modelled by N(e) = No + N\ cos(27re/tj c ). 
The final integral over e can now be done, and at least 
for uj/uj c large compared to No /N't, the result is 



Aa xx oc — sin(27rw/o; c ), 



(2) 



with a positive coefficient. This form, which resembles 
the derivative of the density of states (dN/de)\ e=u , arises 
because the main contribution is from initial states near 
the center of filled broadened Landau levels, which are 
scattered to empty broadened levels, and the available 
phase space is enhanced or diminished for Ax positive or 
negative, depending on the energy change to modulo oj c 
(see Fig.^l. It is clear from experiment that a xy is nearly 
100 times larger than a xx and is not significantly affected 
by the radiation. Therefore, inverting the conductivity 
tensor yields p xx « p xy a xx , which has the period and 
phase of the oscillations observed in experiment. 

While the above treatment is highly oversimplified, it 
indicates that disorder plays a crucial role and may be 
all that is necessary to obtain the radiation-induced os- 
cillations. This suggests that a diagrammatic (Kubo for- 
mula) calculation of the conductivity, including radia- 



tion and disorder but neglecting electron-electron inter- 
actions, will be sufficient to reproduce the effect. This 
calculation is presented below. 

In the Landau gauge and the rotating- wave approxima- 
tion (we neglect both counter-rotating and guiding-center 
terms in the coupling to radiation), the Hamiltonian is 

H = '^2 m J c cl lk C nk + C nk c n'k'V n ^-n',k' 
k,n k.n;k' ,n' 

, e ^ I— ( t —tort i t iut \ fn\ 

+ —^}^^ n { C nk C n-lMe +<_l,fc<^fcC J (3) 
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where n is the Landau level index, k is the y-component 
of momentum, c n k is the electron annihilation operator, 
V n> k;n' ,k' are the matrix elements of the disorder po- 
tential Vj mp (r), and E is the magnitude of the electric 
field component of the microwave radiation. Following 
the calculation of Ando 5] for the zero-radiation case, 
we shall include disorder within the self-consistent Born 
approximation (SCBA) and assume Vimp(r)Vi m p(r') = 
(2<y/m*)8(r — r') where 7 = l/2r is the elastic scat- 
tering rate in zero magnetic field. This assumption of 
(5-corrclatcd disorder, while somewhat inappropriate for 
the long-ranged impurity potentials associated with mod- 
ulation doping, significantly simplifies the problem by 
eliminating the momentum dependence of the self-energy. 
Yet, as we shall see, it manages to capture the impor- 
tant physics rather well. To all orders in the disorder 
and radiation, the Green's functions are given by the di- 
agrams in Fig. |2Ja). Since the radiation terms connect 
neighboring Landau levels, evaluating these diagrams for 
the retarded or advanced Green's functions, G^j^(ti,t 2 ), 
yields a matrix equation in Landau level indices n and 
to. (The Green's functions are proportional to an iden- 
tity matrix in the momentum label k, which is there- 
fore dropped.) The presence of radiation renders this 
an inherently non-equilibrium problem which requires 
the approach of Kadanoff and Baym and Keldysh (see 
Refs. 0, H for details). This involves the "lesser" 
Green's function, G^ m (ti,t 2 ), which is the only place 
where the distribution of electrons enters, and which is 
evaluated from the same set of diagrams via a slightly 
more complex procedure. It is convenient to define 
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inujtx ^<R,A,< 1, . \ —imujt2 

- G nm (t 1 ,t 2 )e 



(4) 



where t = h-t 2 ,T = (h+t 2 )/2, and G^'<(z,T), are 
the Fourier transforms in t. After some manipulation, we 
obtain the following system of equations: 

[z — n(u) c — u>) — Yi R ' A (z + nui) ± i8) 5 np 



■— /=■ (V n Vp+l + VP°n+l,p) 



Gpm( z ) — 8 n m (5) 
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tions, it is straightforward to calculate the conductivity 
from the Green's functions. We find that 



+ 



(b) 

dc 



+ 



dc 



+ 





s dc dc -Sr ^"N* dc 



FIG. 2: Diagrams for (a) Green's function and (b) polariza- 
tion bubble including radiation and disorder within SCBA. 
Disorder lines are dotted. Photon lines are curvy. The ver- 
tex of the polarization bubble is dressed with ladders of dis- 
order lines. When the no-photons-crossed-by-disorder-lines 
conserving approximation is made, diagrams like those in (c) 
are included while those in (d) are neglected. 



G< m (z) =J2 G* p {z) [n F (z+poj)2 t 5+Z<(z+pu;)] G* m (z) 
p 



(6) 



^■A,< {z) = T±J2G^<(z-nu) (7) 

7T * — ' 
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where S — ► + . The Eqs. I r >l7l) have been displayed for 
the steady state case, in which all Green's functions are 
independent of the average time, T. The general T- 
dependent situation leads to more complicated expres- 
sions which can be obtained as discussed in Ref.0. Given 
the Green's functions, the conductivity is obtained via 
the Kubo formula, dij — — limn— >o I m (fi)/fi, where 
nf5 is the retarded current-current correlation function, 
and i, j = {x, y}. Diagrammatically, this means evaluat- 
ing the polarization bubble in Fig. Efb) where the vertex 
is dressed by ladders of impurity lines. 

An exact solution is complicated by the self-consistent 
nature of the above equations. We shall simplify by ne- 
glecting vertex corrections to the polarization bubble. To 
do so within a conserving approximation, we must also 
neglect all diagrams in which impurity lines cross photon 
insertions (see Fig.[2Ic,d)). Note that this is equivalent to 
replacing the fully self-consistent E's in the above equa- 
tions with E's calculated in the absence of radiation but 
still self-consistent in the disorder. Thus, if we define 
E(z) to be the self-energy calculated by Ando for the 
zero-radiation case, then Y, R (z) — ► E(z), E A (z) — > E*(z), 
andE < (z) — > — 2inp{z)hrxS{z). This approximation also 
has the virtue of mimicking energy loss mechanisms so 
that a steady state T-independent solution does indeed 
exist 0, and can be obtained directly from Eqs. JHEJ- 
The radiation enters only when the Green's functions are 
calculated from the E's. In the absence of vertex correc- 
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^2 y/(n+l)(m+l) J dz 



(8) 



and obtain the resistivity by inverting the cr-matrix. 

Within the approximations discussed above, and us- 
ing parameter values appropriate to the experiments of 
Ref.0, we have calculated p xx numerically to all orders in 
radiation and disorder. Results are presented in Fig. 
where we plot p xx versus l/u> c for fixed w. We consider 
three values of radiation intensity (power per unit area 
in units of m*u> 3 ): 1 = (dark), I = 0.0034 (moderate 
intensity), and I = 0.0115 (high intensity). The dark re- 
sistivity exhibits only the familiar SdH oscillations which 
have period l//z and decay away as lu c becomes small 
compared with the temperature. For moderate intensity, 
we find a pronounced radiation-induced oscillation of p xx . 
In agreement with experiment, the period of oscillation is 
1/uj and minima are found nearuj/uj c = integer-)- 1/4. We 
note, however, that the 1/4 phase shift is not universal, 
varying between and 1/2 depending upon disorder and 
intensity. A more robust feature is the presence of zeros 
of the oscillation at the integer values of u)/ui c . Unlike the 
SdH oscillations, the radiation- induced oscillations are 
not sensitive to the temperature-broadening of the elec- 
tron distribution about the Fermi level. T-dependence 
therefore derives from the the T-dependence of the scat- 
tering, a feature absent from the present calculation due 
to our neglect of interaction effects. The magnitude of the 
radiation-induced oscillations can easily exceed the dark 
resistivity, leading to regions of negative total resistivity. 
This is reasonable for a non-equilibrium system, and a 
similar effect has been observed experimentally in semi- 
conductor superlattices 0, but it seems that some ad- 
ditional physics is required for these negative-resistivity 
minima to become the zero-resistance states observed in 
the present experiments. For high intensity radiation, ad- 
ditional features appear which have yet to be discovered 
experimentally. These likely correspond to multi-photon 
processes by which an electron absorbs m photons and is 
promoted by n Landau levels. 

Since the simplifying approximations discussed above, 
employed to make the calculation tractable, are not well 
controlled, we cannot expect our results to be quantita- 
tively accurate. However, the fact that even this sim- 
plified calculation captures the qualitative form of the 
experimental results encourages us that we have identi- 
fied the essential physics of the problem. Due to the as- 
sumption of (^-correlated disorder, even the zero-radiation 
(dark) resistivity, p xx , does not agree quantitatively with 
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FIG. 3: Calculated radiation-induced resistivity oscillations. 
We plot p xx vs 1/Uc at fixed ui for [i — 50a>, fcsT = w/4, 
7 = 0.08a;, and three values of radiation intensity (in units 
of m*u 3 ): 1 = (dark), I = 0.0034 (upper panel), and I = 
0.0115 (lower panel). For computational purposes, the energy 
spectrum is cutoff at 20 Landau levels above and below the 
chemical potential. The high-frequency oscillations seen at 
small 1/wc are the familiar SdH oscillations with period l//u. 



experiment. As the calculation neglects the difference be- 
tween the transport lifetime, r tr , and the single particle 
lifetime, r, it overestimates p xx by a factor of r^/r, equal 
to about 50 in the present case. Since the dark Hall re- 
sistivity, p xy , is approximately independent of disorder 
in the -B-field regime of interest, the calculated value of 
p x is correct, but effectively too small compared to p xx . 
This quantitative error in the dark resistivity has two im- 
portant consequences. When we calculate the radiation- 
induced change in the Hall resistivity, Ap xy , we find that 
it is also oscillatory and of comparable magnitude to 
Ap xx . Since p xy is effectively too small, we find a total 
Hall resistivity that exhibits noticeable radiation-induced 
oscillations. For analogous reasons, while our calcula- 
tions yield oscillations in a xx which grow linearly with ra- 
diation intensity, the intensity-dependence is muted when 
we invert the conductivity matrix to obtain p xx . How- 
ever, correcting for the missing factor of 50, it is clear 
that the effect of radiation on p xy is negligible and that 
the oscillations in p xx grow linearly with intensity, which 
is what was observed experimentally. 



In conclusion, we understand the radiation-induced 
magnetoresistance oscillations recently observed in the 
2DEG to be a consequence of photo-excited disorder- 
scattered electrons contributing to the dc resistivity in 
an oscillatory manner inherited from the broadened Lan- 
dau level structure of the energy spectrum. Our dia- 
grammatic calculation yields a radiation-induced oscilla- 
tion of p xx that has the correct period and phase. We 
find that for radiation of sufficient magnitude, the am- 
plitude of the oscillation can be made large enough that 
the total resistivity is negative in the vicinity of the min- 
ima. This effect is reminiscent of the experimentally 
observed zero-resistance states but indicates that addi- 
tional physics, perhaps localization or electron-electron 
interactions, is required to understand this aspect of the 
observed phenomenon. Our calculations at high inten- 
sity reveal multi-photon effects, the observation of which 
should be a goal for future experiments. More advanced 
calculations, going beyond the simplifying approxima- 
tions employed herein and including the interaction ef- 
fects responsible for energy relaxation and temperature 
dependence, are left for future work. 
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After submitting this Letter for publication, we learned 
that a related physical picture and a calculation of the 
conductivity to second order in the radiation, with a pre- 
diction of negative conductivity, were previously reported 
by Ryzhii et al. 
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